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Abstract 

The SELEX Collaboration has reported a very large isospin splitting of doubly charmed 
baryons. We show that this effect would imply that the doubly charmed baryons are very 
compact. One intriguing possibility is that such baryons have a linear geometry Q — q — Q 
where the light quark q oscillates between the two heavy quarks Q, analogous to a linear 
molecule such as carbon dioxide. However, using conventional arguments, the size of a 
heavy-light hadron is expected to be around 0.5 fm, much larger than the size needed to 
explain the observed large isospin splitting. Assuming the distance between two heavy 
quarks is much smaller than that between the light quark and a heavy one, the doubly 
heavy baryons are related to the heavy mesons via heavy quark-diquark symmetry. Based 
on this symmetry, we predict the isospin splittings for doubly heavy baryons including 
Sec, Sbf, and She. The prediction for the Sec is much smaller than the SELEX value. On 
the other hand, the Sf,h baryons are predicted to have an isospin splitting as large as 
(6.3 ± 1.7) MeV. An experimental study of doubly bottomed baryons is therefore very 
important to better understand the structure of baryons with heavy quarks. 
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1 Introduction 



A key prediction of QCD is the existence of baryons with two or three charm or bottom quarks. 
Several years ago, evidence for the hadroproduction of five different baryons with two charm 
quarks was reported by the SELEX Collaboration at Fermilab [H El [3l HI [5l [6] . Two singly 
charged states (3443) and S+ (3520) were observed in the A+i^~7r"'" mass distribution, and 
three doubly charged states 5++ (3460), H++(3541) and 5++ (3780) were observed decaying 
into A+K-7r+7r+ PEIIS]. The H+(3520) was also observed in the pD+K~ [4J and H+vr+vr^ 
final states [6]. An analysis of the helicity angular distribution support the assignments that 
the (3443) and H^"^(3460) form an isospin doublet, and the (3520) and H^''"(3541) form 
another. The preliminary isospin mass splittings were reported to be 17 MeV and 21 MeV, 
respectively [3]. This observation is very puzzling because such values are much larger than 
all isospin mass splitting of hadrons known so far. For instance, the mass difference between 
the proton and neutron is — rup = 1.29 MeV, and between the charged and neutral D 
mesons is Mjj+ — M^o = 4.77 ± 0.10 MeV [7]. The largest isospin splitting ever observed is 
the double-strange baryons M=- - M=o = 6.85 ± 0.21 MeV [7]. 

Recently, the mass of the lowest H^"*" state was updated from 3460 to 3452 MeV [8]. 
Although the isospin splitting is decreased from 17 MeV to 9 MeV for the lower doubly 
charmed baryon isospin-doublet, it is still larger than all the other known isospin splittings. 
It is thus interesting to see whether it is possible to obtain the observed rather large values 
from known physical principles with controlled uncertainty. In Section [21 it will be shown that 
the SELEX observations would imply the Hcc to be very compact, which, however, cannot be 
understood by any known mechanism of the strong interactions. 

Predictions for the isospin splittings of doubly heavy baryons will be presented in Sec- 
tion [3] based on the conventional assumption that the two heavy quarks constitute a compact 
diquark. We then can apply the ansatz of heavy quark-diquark symmetry. Our predictions 
for the doubly charm baryons are similar in magnitude to the isospin splittings for other 
hadron isospin multiplets, but considerably smaller than the SELEX data. We will also give 
predictions for the isospin splittings of the ccccq and bbbbq pentaquark states based on a heavy 
quark- "quadra-quark" symmetry in Section HI 

SELEX used pi^h = 600 GeV/c 7r~,S~, and proton beams on a nuclear target to pro- 
duce the doubly charm baryons. A striking feature of the SELEX measurements is the fact 
that the observed doubly charmed baryons are all produced at xp > 0.1 (SELEX only has 
sensitivity in that region); i.e., at a significant fraction of the projectile momentum. This 
is consistent with the ISR measurements of the Ac [9j and the A;, [IDJ at high xp, as well 
as NA3 measurements at CERN [11] which showed that two J/ip^s are hadro-produced at 
high Xi? in pion-nucleus collisions; in fact, each ttA — t- J/ijjJ/ipX event measured by NA3 
has four charmed quarks with a flat longitudinal momentum distribution for xp > 0.4. [12j. 
The SELEX and NA3 measurements cannot be explained if the heavy quarks only arise from 
gluon splitting; however, this is a natural consequence of the existence of intrinsic heavy 
quarks in the projectile [131 [T4l [151 116j . such as the rare \uudcccc) Fock state in the proton 
or the \udcccc) Fock state in the vr". Since the momentum distribution in such Fock states is 
maximized at low invariant mass, all of the quarks tend to have the same rapidity and small 
transverse momentum. The heavy quarks have the maximum momentum fractions in such 
configurations since equal rapidity implies Xi oc \J'mf + k'j_^. The doubly charmed ccq baryons 
are then formed in a collision by the coalescence |17[ [T8] of the comoving heavy quarks with 
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a light quark of the projectile — the domain where the wave function of the produced doubly 
charm hadron is maximal. This mechanism also explains why doubly charmed baryons are 
not readily produced in e+e~ annihilation; in that case it is rare for the two charmed quarks to 
be in the same kinematic domain. The intrinsic charm mechanism also accounts for the non- 
factorized nuclear-target dependence [TH1[2D] of J/ifj hadroproduction [TTIIIB]. It also points 
to the high xp domain of hadroproduction as the best kinematic region to search for heavy 
hadron systems in general. Thus the best opportunity to create superheavy hadrons and test 
their properties is in hadron-hadron collisions at high x p using the intrinsic heavy quark Fock 
state mechanism — for example, at the LHCb, or at future fixed-target experiments using 
the 7 TeV LHC beam. 

The last section contains a brief summary of our results. 



2 Implication of a large isospin splitting 

It is instructive to ask the question: What does a large isospin splitting imply for the doubly 
heavy baryons? Isospin splittings originate from two sources — the u and d quark mass dif- 
ference as well as electromagnetic contributions. The interference pattern of the two different 
contributions to the mass differences can be easily understood. The repulsive (attractive) 
Coulomb interaction gives positive (negative) contribution to the electromagnetic (em) self- 
energy of the baryons, so that the baryon with more absolute electric charge has more em 
self-energy. The sign of the quark mass difference contribution reflects the fact that the down 
quark is heavier than the up quark. Hence, for the and Hfec, the interference is destructive 
while for the Hf,;, it is constructive. The sign of the em contribution to M-++ — is the 
same as the SELEX data, however, the em and the quark mass term interfere destructively 
substantially reducing the em effect. Thus, in order to quantitatively understand the SELEX 
data an unusually large em contribution to the mass differences is necessary. 

In the case of heavy particles the em effect is Coulombic since the magnetic contribution 
is negligible. To quantify the em contribution, one may employ the Cottingham formula to 
analyze the contribution of virtual photons. It can be used to relate the em self-energy to the 
em form factor of a particle, see e.g. [21j. Then the em self-energy is given by |21] (neglecting 
the inelastic contributions) 

aQ^ fd^l^r. , 2m2 



where a = 1/137.06 is the fine structure constant, Q is the total electric charge in units of the 
proton charge, and Gpit) is the Fourier transform of the charge distribution of the particle. 
Taking a dipole distribution with a mass parameter m in units of GeV (which is sufficient for 
our purpose) 

^^(*) = (l-t/m^y^ (2) 
one can perform the integration analytically, and gets 

= i"^'"^- (3) 

For a first estimate, let us consider the ccq as a two-particle system with charge 4/3 and 
Cq. Therefore, the em contribution to the isospin splitting of the doubly charmed baryons 
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6" = M^++ - with AQ2 = 3 is 

^cc ^cc 

5|P = —am = 0.0034 m. (4) 
32 

The mean square radius of a particle is obtained by taking the first derivative of its em form 
factor. For a heavy particle with the em form factor given by Eq. ([2]), it is 



dt 



t=o m 



12 

-• (5) 



Since the and contain ecu and ccd quarks, respectively, the contribution to 5-^^^ 
from the quark mass difference is negative, i.e. 5=™°^ < because the u quark is lighter than 
the d quark. Hence, the em contribution must be larger than the total isospin splitting of 
the doubly charmed baryons, > fe^^. If we take 9 MeV as the isospin splitting, one gets 
m > 2.65 GeV from Eq. dU, and 



2) < 0.26 fm. (6) 

This value is much smaller than the typical size of a hadron containing light quark(s). In 
fact, it is of similar size as the distance between the two heavy quarks ~ l/(mQf), with mq 
and V the mass and the velocity of the heavy quark, respectively. If we use a larger isospin 
splitting, e.g. 17 MeV, instead, the resulting (r^) < 0.14 fm is even smaller. 

Therefore, we conclude that a large isospin splitting would imply the doubly heavy baryon 
to be very compact — the larger the splitting, the smaller the size. This conclusion can be 
easily understood because a large isospin splitting would mean a large em self-energy which, 
being proportional to (1/r), in turn would mean a small size of the doubly heavy baryon. 

One intriguing possibility is that doubly charm baryon states have a linear geometry 
Q — q — Q where the light quark q oscillates between the two heavy quarks Q, analogous 
to a linear molecule such as carbon dioxide CO2 = O — C — O. In this case the overall 
size of the baryon would be relatively small. A lattice gauge theory investigation of this 
possibility would be interesting. However, we are not aware of a mechanism in quantum 
chromodynamics (QCD) which can keep the light quark in line and close to the heavy quarks. 
In particular, if we take a Coulomb plus linear potential as the interquark interaction, the 
distance between the light quark and a heavy quark scales as (cmcons)""'^^'^ ~ 1/Aqcd with 
the flavor-independent string tension a being the coefficient in front of the linear potential 
and w-cons the constituent light quark mass. The numerical value for ^/o^ is about 430 MeV 
from Regge trajectories of light mesons and also from heavy quarkonia spectrum, see e.g. [22j . 
Thus the size of a heavy-light hadron is expected to be around 0.5 fm, much larger than the 
size needed to explain the observed large isospin splitting. 



3 Isospin splittings of doubly heavy baryons — a quantitative 
analysis 

In the following we will assume that the distance between the light quark and a heavy quark is 
much larger than the distance between the two heavy quarks inside the doubly heavy baryon. 
Based on this conventional assumption, it was proposed that there is a heavy quark-diquark 
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symmetry [23] which relates a doubly heavy baryon containing two heavy quarks to a heavy 
meson containing a heavy anti-quark. 

The distance between the two heavy quarks in a doubly heavy baryon is characterized 
by TQQ ~ l/{mQv), with v being the heavy quark velocity. It is much larger than the 
length scale for the light quark r^g ~ I/Aqcdj as already discussed in the previous section. 
Hence, one may perform an expansion in rgg/rgQ with a controlled uncertainty. In the heavy 
quark limit, only the leading term is necessary, which means the diquark formed by the two 
heavy quarks is point-like. Furthermore, the diquark has the same color charge as a heavy 
anti-quark. Hence, there is an approximate U(5) symmetry relating the spin-1/2 heavy anti- 
quark and the spin-1 heavy diquark, called heavy quark-diquark symmetry [23]. Using this 
symmetry, the doubly heavy baryons can be studied by relating them to the heavy mesons. 
In this section, we will calculate the em as well as quark mass difference contribution to the 
isospin splitting of doubly heavy baryon masses at next-to-leading order (NLO) in the chiral 
expansion. This is the lowest order at which isospin breaking operators appear. In view of 
the present knowledge of the masses of the doubly heavy baryons, this should suffice. Both 
effects can be taken into account systematically up to a given order using chiral perturbation 
theory with virtual photons \24:\ [25] . We will use the formalism proposed in Ref . [26] , which 
combines heavy mesons and doubly heavy baryons into the same field, and construct the NLO 
chiral Lagrangian which is responsible for the isospin mass splittings. 

The super-flavor multiplet of heavy mesons and doubly heavy baryons can be collected 
into a single 5x2 matrix field I26j, which can be written in components as 

^a,/i/3 = -f^a,a/3 + Ta,i/3 (7) 

where a = u,d,s is the light flavor index, fi runs from 1 to 5, a, /3 = 1, 2 and i = 1, 2, 3. The 
fields for the heavy mesons. Ha, and and doubly heavy baryons, Ta, are given by 

Ha,al3 = Pa ' ^al3 + Pa^ajS, 
Ta,ip = \/2 [r^a^ip + -j=r.a,a(yl,l3^ , (8) 

where Pa and Sa are the fields annihilating the vector (pseudoscalar) heavy mesons and 
the spin-1/2 (3/2) doubly heavy baryons, respectively. The field for the spin-3/2 baryon is 
constrained by H* j^o"^^ = 0. 

One can construct the effective Lagrangian for the mass terms assuming the heavy quark- 
diquark symmetry. At leading order, there is no isospin splitting within the same multiplet 
as can been from the Lagrangian constructed in Ref. |26) . At NLO, the Lagrangian relevant 
for the isospin and SU(3) mass differences reads 

Asv = -cTr [nl'Hi,{x+)ba - doTi niq'Hh{Q+)ba 

-TV {niH, [diiQl - Q-)ba + d2{Q+{Q+))ba] } , (9) 

where the heavy-flavor charge operator q is defined as qHa = qqHa and qTa = —2qQTa, with 
qq being the charge of the heavy quark in a heavy meson. It is similar to the isospin breaking 
terms in the Lagrangians constructed for heavy mesons [27] and singly heavy baryons [28j. 
The operators x+ and Q± contain the Goldstone boson fields which are needed for higher 
order calculations 



uhn^ + uxn, Q^ = Uu'Qu±uQu^), (10) 
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where u = VU, U = exp ^-v/2z0/F^^, with the pion decay constant, and (f) cohects the 
pseudoscalar mesons, 



K- K° 



(11) 



The Ught quark mass and charge matrices are 

X = 2Bo- diag {ruu, rud, rUs} , Q = e • diag . (12) 

The Ught quark mass and em contributions to the heavy meson and doubly heavy baryon 
masses can be easily worked out using this Lagrangian. At NLO, the Goldstone bosons are 
not needed since the chiral loop corrections to the hadron masses start from higher order. In 
this case, physically, the di term describes the effect of the virtual photons coupled to the 
light quark. Because there is only one light quark in both heavy mesons and doubly heavy 
baryons, these virtual photons merely contribute to the self-energy of the light quark. So 
they can be absorbed into a redefinition of the quark masses 

m„ = m. + ^, -^W=-<^W + ^- (13) 

There is no d2 term because ((5+) = {Q) vanishes with the charge matrix given in Eq. (jl2p . 
We remark that we are well aware of the subtleties concerning em corrections to quark masses 
(see e.g. Ref. [29] ) but these can be ignored to the accuracy we are working. 

Therefore, at NLO there are only two parameters describing the mass corrections of the 
heavy mesons and doubly heavy baryons. One is for the light quark mass difference, and 
the other one is for the em effects which originate from the virtual photons exchanged be- 
tween the light quark and the heavy quarks. These two effects can also be parameterized 
by two parameters in quark models once neglecting the spin-dependent interactions which 
are suppressed by l/mq, see e.g. [30l|3l]. The explicit expressions for the mass corrections 
can be found in the Appendix. The former parameter can be determined through the mea- 
sured SU(3) mass splitting between the heavy mesons, and the latter one can be determined 
from the isospin mass splitting between the charged and neutral heavy mesons [27]. Defin- 
ing c = 4c-Bo("^s — "^d) and d = d^e^ /3, using the mass differences among the pseudoscalar 
charmed mesons D^,D~^,Df [7], one gets 

c = (98.99 ± 0.30) MeV. (14) 

The value of d can be extracted via 

d = -[{M^+ - Md+)X - {Md+ - M^o)] = (-1.05 ± 0.16) MeV, (15) 
2 * 

where A = (md — rhu)/{ms — fhd) = 0.027 it 0.003 is calculated from the recent FLAG 
average [32^j. This value for A is consistent with the latest precise lattice determinations of 
the light quark masses, A = 0.029 it 0.002 [33]. Note that in both cases the value for A refers 
to the masses without the em shift. In what follows, we will ignore this difference, as it is 
expected to be a minor effect. Therefore, these values correspond to 

(M^+ - Mbo)''*™^s = (2.67 ± 0.30) MeV, (M^+ - M^o)^™ = (2.10 ± 0.32) MeV, (16) 
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M„+ -M=o 


EM 


4.2 ±2.3 


4.0 ± 1.5 


1.6 ± 1.1 


Strong 


-2.7 ± 1.5 


2.3 ±0.8 


-2.5 ± 1.4 


Total 


1.5 ±2.7 


6.3 ± 1.7 


-0.9 ± 1.8 


Ref. [31J 


2.3 ± 1.7 


5.3 ± 1.1 


-1.5 ±0.9 



Table 1: Predicted isospin splitting of the doubly heavy baryons. The electromagnetic and 
quark mass difference contributions are given in the second and third rows, respectively. The 
final results are shown in the fourth row. The results in Ref. [31] are given in the last row for 
comparison. All values are given in units of MeV. 



where the superscript "strong" and "em" denote the effects from the u, d quark mass difference 
and the virtual photons, respectively. If we use the mass differences among the pseudoscalar 
bottom mesons, M^o - Mb = 87.0 ± 0.6 MeV and M^o - Mb+ = 0.33 ± 0.06 MeV 0, then 
we get different values 

c = (87.0 ± 0.6) MeV, d = (-2.0 ± 0.3) MeV, (17) 

which corresponds to the decomposition 

(M^o - MB+f^'°'^^ = (2.3 ± 0.3) MeV, (M^o - Mb+)^™ = -(2.0 ± 0.3) MeV. (18) 

The values agree within uncertainties — the differences between these two sets of values may 
be understood as flavor symmetry breaking corrections of order ©(Aqcd/^t-c)- 

Using these parameter values, the mass difference between the H++ and can be easily 
obtained. One can also get the isospin splittings for the doubly bottom baryons Stb and 
the bottom-charm baryons 5^^. All of the predictions are listed in Table [TJ To minimize 
the uncertainty from heavy quark flavor symmetry, the results for the Hcc and H^f, are given 
using the values of c and d extracted from the charmed and bottom mesons, respectively. 
The results for the Hfec cover the values obtained using both sets of parameter values. Now 
let us discuss the other possible uncertainties. Except for the ones from c and d, there are 
still uncertainties from neglecting higher order counterterms and loops in the chiral expansion. 
Since our predictions concentrate on the isospin splittings, the relevant corrections from higher 
order terms in the chiral expansion should be of order ©(Mtt/A^) ~ 15%, with M^- the pion 
mass and A^ ~ 1 GeV the chiral symmetry breaking scale. In addition, there should also 
be corrections to the heavy quark-diquark U(5) symmetry. These corrections should be of 
order 0{rQQ/rqQ) = O {Aqqd / {"rnQv)) which describes the relative size of the neglected heavy 
quark separation with respect to the distance between the light quark and the heavy diquark. 
Conservatively, we take O {Aqcd / {mQv)) ~ 50%, 30% and 50% for the Ecc, '^bb and Ebc 
baryons, respectively. The same isospin splittings were also calculated in an approach based 
on a parameterization inspired by heavy quark effective theory and utilizing some data to fix 
the parameters [31]. For comparison, their results are given in the last row. 

One can get both spin-3/2 and 1/2 doubly heavy baryons from binding a spin-1 heavy 
diquark and a light quark. Because the spin of the diquark is the same in both cases, they are 
related to each other by the heavy quark spin symmetry, and have the same isospin splittings 
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as given in Table [H Corrections to the spin symmetry are suppressed by Aqcd / itt-q ■ These 
corrections are expected to be smah as confirmed by a comparison with the results of a quark 
model [3D] which takes into account the spin-dependent interactions, namely M^- — M^^o = 

^bb ^bb 

6.24 ± 0.21 MeV or 6.4 ± 1.6 MeV using different inputs, which is quite close to ours. 

4 Isospin splittings of quadratically heavy pentaquarks 

The analysis can be extended to pentaquarks containing four heavy quarks and one light 
quark. One should notice that the size of the four heavy quark cluster is not four times larger 
than the distance between two heavy quarks, tqq. Being in an 5- wave, if the four quarks are 
of the same flavor, they should be spatially symmetric. Hence, the size of the cluster is the 
same as rgg ~ l/('mQv), which is again much smaller than 1/Aqcd in the heavy quark limit. 
Moreover, in a pentaquark, the four heavy quarks are in a fundamental representation of the 
SU(3) color symmetry group, which is the same as one quark. Hence, to a first approximation 
they can be treated as one object, to be called quadra-quark in the following. For the four 
heavy quarks of the same flavor being in an 5- wave, Fermi statistics requires their spin wave 
function to be symmetric. Hence, the quadra-quark is a spin-2 object. Analogous to the U(5) 
symmetry between the heavy diquark and heavy anti-quark, the symmetry for the heavy 
quadra-quark and heavy quark is U(7). One may find an interesting phenomenology for the 
pentaquarks using the U(7) symmetry. Here, we are only interested in the isospin splittings. 
It is easy to find the following relations to lowest order in isospin breaking: 

Kcccd - Mccccu = 4(M^+ - MdoT"^ + (M^+ - M^o)^*'-°°g = (11 ± 5) MeV, 

M,,,,j - Mmbu = 4(Mbo - Ms+ + (M^o - = - (6 ± 3) MeV. (19) 

The splittings are large, but certainly more of an intellectual curiosity at present. 

5 Summary 

Using the Cottingham formula to compute the Coulomb electromagnetic shift, we have shown 
that the large SELEX value of the isospin splitting of the states implies that double charm 
baryons are very compact; i.e., the light quark must be very close to the two heavy quarks. 
A novel possibility is that the quarks in the doubly charm baryons are arranged as a compact 
state c—g — c with a linear geometry. This possibility could be checked using lattice gauge the- 
ory simulations. However, the infrared behavior of the light quark is expected to be governed 
by the non-perturbative confining interaction, and thus the size of any hadron containing a 
light quark should be of order 0(1/Aqcd)- A conventional approach exploiting this is based 
on quark-diquark symmetry. It allows us to predict the isospin splitting for doubly heavy 
baryons Hcc, '^bb and She at NLO in the chiral expansion. These predictions for the doubly 
charm baryons give isospin separations much smaller than the SELEX measurements. There- 
fore, the compactness implied by the SELEX data appears to call for a significant violation 
of heavy quark-diquark symmetry — today no mechanism is known that can provide this. 
However, it should be noticed that among all the four Hcc states in the two reported isospin 
doublets only the mass of the H+ (3520) has been measured with certainty. 

In order to resolve the discrepancy between the experiment and theory, further precise 
experimental investigations are clearly needed. Our prediction for the isospin splitting is 
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6.3 lb 1.7 MeV based on heavy quark diquark symmetry. It is of similar size as M^- — M=o 
which contain two strange quarks. (Because nig < Aqcd; we have refrained from making 
use of the same method to calculate the isospin splittings containing strange quark(s).) Any 
configuration that leads to isospin splittings as large as those reported by SELEX for the 
doubly charmed baryons would lead to significantly larger splittings for the doubly bottomed 
baryons, because the em self-energy would be much larger. 

We have also made predictions for the isospin splittings of the pentaquarks ccccq and bbbbq 
and have found that the value for the ccccq is as large as (11 it 5) MeV using a generalization 
of heavy quark diquark symmetry. 

As we have noted, the best chance to create these super-heavy hadrons and test their 
properties is in hadron-hadron collisions at high xp using the intrinsic heavy quark Fock 
state mechanism — for example, at the LHCb, or at future fixed-target experiments using 
the 7 TeV LHC beam. 
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A Expressions for the mass corrections at NLO 

In the appendix, we give explicit expressions for the NLO corrections of the masses of the 
heavy mesons and doubly heavy baryons. In the heavy quark limit, the QCD Lagrangian 
does not depend on the heavy quark mass, which results in the heavy quark flavor symmetry, 
see e.g. |34) . In the following, the formula are given in terms of general heavy quark flavors, 

2 2 

Qu : AcBorhu + -eqdoe'^, Q1Q2U : 4c5om„ - -(eg, + egj^oe^, 

Qd : AcBofhd - ^eqdoe'^, QiQ2d : 'icB^rhd + ^(eg^ + egjdoe^, 

Qs : AcBorhs - ^cqdoc^ , Q1Q2S : icBorhs + ^{cq^ + eQ^)dQe^ , (A.l) 

where Q(i,2) represents the heavy quark flavor, and eg^^ its charge in unit of the elementary 
charge e with e > 0. 
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